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9^' ABSTRACT 

<: 

CSj ' The cosmological magnetic field is one of the important physical quantities 

^ ' which affect strongly the cosmic microwave background (CMB) power spectrum. 

t> ' Recent CMB observations have been extended to higher multipoles I >1000, and 

^3 ■ they resultantly exhibit an excess power than the standard model prediction in 

cosmological theory which best fits the Wilkinson Microwave Anisotropy Probe 
(WMAP) data at lower multipoles I <900. We calculate the CMB temperature 
anisotropies generated by the power-law magnetic field at the last scattering 
surface (LSS) in order to remove the tension between theory and observation at 
Oh 1 higher multipoles and also place an upper limit on primordial magnetic field. In 

O ' our present calculation we take account of the effect of ionization ratio exactly 

without approximation. This effect is very crucial to precisely estimate the effect 
of the magnetic field on CMB power spectrum. We consider both effects of 
the scalar and vector modes of magnetic field on the CMB anisotropies, where 

^ ■ current data are known to be insensitive to the tensor mode which we ignore 

in the present study. In order to constrain the primordial magnetic field, we 
evaluate likelihood function of the WMAP data in a wide range of parameters of 
the magnetic field strength \~B\\ and the power-law spectral index rig, along with 
six cosmological parameters in flat Universe models, using the technique of the 
Markov Chain Monte Carlo(MCMC) method. We find that the upper limit at 2a 
C.L. turns out to be |B>| < 3.9nG at 1 Mpc for any n# values, which is obtained 
by comparing the calculated result including the Sunyaev-Zeldovich(SZ) effect 
with recent WMAP data of the CMB anisotropies. 
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1. Introduction 

Temperature and polarization anisotropies in CMB provide very precise information 
of the physical processes in the early Universe. However, recent new CMB data sets from 
the WMAP (Bennett, et al. 2003), the Arcminute Cosmology Bolometer Array Receiver 
(ACBAR; Kuo et al. 2004), and the Cosmic Background Imager (CBI; Mason et al. 2003) 
have indicated a potential discrepancy between theory and observation at higher multipoles 
I > 900. The best-fit cosmological model to the WMAP data predicts the power spec- 
trum which shows appreciable departure from those observed in balloon and interferometer 
experiments. This discrepancy cannot be explained by taking account of tuned standard 
cosmological parameters. 

One possible interpretation of an excess power at high multipoles is a manifestation 
of the re-scattering of CMB photons by hot electrons in clusters known as the SZ effect 
(Sunyaev and Zeldovich 1980). Although there can be no doubt that some contribution from 
the SZ effect exists in the observed angular power spectrum, it has not yet been established 
conclusively that this is the only possible interpretation (Aghanim et al. 2001) of the small 
scale power. Indeed, the best value of the matter fluctuation amplitude to fit the excess 
power at high multipoles is near the upper end of the range of the values deduced by the 
other independent method (Bond et al., 2002; Komatsu and Seljak 2002). Toward the 
solution of this bothersome problem, it has recently been reported that the bump feature in 
the primordial spectrum gives the better explanation for both CMB and the matter power 
spectra at small scales (Mathews et al., 2004). 

The inhomogeneous cosmological magnetic field generated before CMB last scattering 
epoch is also a possible candidate to reconcile the tension in such higher multipoles. It ex- 
cites an Alfven-wave mode in the baryon-photon plasma in the early Universe and induces 
small rotational velocity perturbations. Since the mode can survive on scales below those of 
the Silk damping during recombination (Jedamzik et al., 1998; Subramanian et al., 1998), it 
could be a new source of the CMB anisotropies in such small scales. Analytic expressions of 
temperature and polarization angular power spectra in rather larger angular scales (I < 500) 
based upon the thin LSS approximation were derived for both vector and tensor modes 
(Mack et al., 2002). Subramanian and Barrow (2002) considered the vector perturbations 
in the opposite limit of smaller angular scales. Too strong magnetic field in the early Uni- 
verse, however, has a possibility of conflicting with the cosmological observations currently 
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available. The combination of those studies and current observations places a rough bound 
on the strength of primordial magnetic field |B A | < (1.0 — 10)nG. 

In order to compare the theoretical CMB anisotropics calculated by including a primor- 
dial magnetic field, with observations much more precisely, we need to perform numerical 
calculations of the fully linearized equations with the use of realistic recombination history 
of the Universe. In particular, we first need to develop a numerical method to predict the 
theoretical spectrum for intermediate angular scales, in which the analytic approximation 
becomes inappropriate. The numerical study of only the scalar mode was introduced by Koh 
and Lee (2000), and they discussed the sensitivity of CMB to the primordial magnetic field 
strength. Lewis (2004) has recently shown the characteristic of a primordial magnetic field 
effect on CMB in vector and tensor modes, numerically. However, there is no systematic 
study to put constraint on parameters of primordial magnetic field by taking account of both 
scalar and vector modes simultaneously. 

The purpose of this Letter is to explore more completely the effects of a primordial 
magnetic field on the CMB anisotropies, and to place a new limit on the field strength 
together with many other cosmological parameters using the MCMC method in the analysis 
numerically. The present calculation differs from all previous discussions in the following 
important ways. First, we take account of the effect of ionization ratio exactly without 
approximation. Second, we consider the effect of scalar mode of magnetic field more precisely 
than the analysis of Koh and Lee (2000) by introducing the cutoff wave number physically 
consistent with the magneto-hydrodynamic (MHD) study of damping Alfven wave (Jedamzik 
et al., 1998; Mack et al., 2002). The effect of scalar mode of primordial magnetic field changes 
the CMB power spectrum significantly at larger scales for I < 900 (Koh and Lee, 2000), while 
the effect of vector mode, which we also consider in the present study, dominates at smaller 
scales. These effects are very crucial to quantitatively estimate the effects of magnetic field 
on the CMB power spectrum over a wide range of multipoles I for various scales. Third, it 
is for the first time to discuss the primordial magnetic field as a possible solution to remove 
the tension between theory and observation which exhibits in smaller angular scales of the 
CMB anisotropies. 



2. Primordial Stochastic Magnetic Field and Baryon-Photon Fluid 

Before recombination, Thomson scattering between photons and electrons and Coulomb 
interactions between electrons and baryons were sufficiently rapid that the photon-baryon 
system behaves as a single tightly coupled fluid. Since the trajectory of plasma particles is 
bent by Lorentz force in magnetic field, photons are indirectly influenced by the magnetic 
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field through Thomson scattering. 

Let us consider the primordial magnetic field created at some moment during the 
radiation-dominated epoch. The energy density of the magnetic field is treated as a first 
order perturbation in a flat Friedmann-Robertson- Walker (FRW) background cosmology. 
Within the linear approximation, the magnetic field evolves as a stiff source, and therefore 
we can discard all back reactions from the MHD fluid onto the field itself. Also we assume 
that the conductivity of the primordial plasma in the early Universe is infinite before decou- 
pling, which is a very good approximation at the epochs which we are interested in. On this 
assumption we can decouple the time evolution of the magnetic field from its spatial depen- 
dence B(?7, x) = Bo(x)/a 2 for very large scale, where 77 is the conformal time, x is the spatial 
coordinate, and a is the scale factor. We also assume that a background primordial magnetic 
field B is statistically homogeneous, isotropic and random. For such a class of magnetic 
field, the power spectrum can be taken as a power- law P(k) oc A;™ (Mack et al. 2002), where 
k is the wave number in the Fourier space, and ub is the power-law spectral index of the 
primordial magnetic field which can be either negative or positive. 

Evaluating the two-point correlation function of the electromagnetic stress-energy tensor 
of the vector mode, we can obtain the isotropic spectrum 
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for k < k D to a very good approximation (Mack et al. 2002) 1 , where the vector mode 
Lorentz force, L^(k), is given by L^(k) = fcll^^k), B\ = |B A | is the magnetic comoving 
mean-field amplitude obtained by smoothing over a Gaussian sphere of comoving radius A, 
and k\ = 2-n/X (A = IMpc in this paper). In this equation ko is the cutoff wave number in 
the magnetic power spectrum defined by 

2 ng+3 

B\ \ n s+ 5 1 k\ \ "s+ s 



where h is the Hubble parameter in units of lOOkm/s/Mpc (Jedamzik et al. 1998; Mack 
et al. 2002). Since the magnetic field source term IT^^k) depends on the magnetic field 
quadratically, the explicit time dependence of the magnetic stress is given by n«(?7,k) = 
n^)(k)/a 4 . Also the electromagnetic stress-energy tensor of the scalar mode is LT^ ~ 



1 This has similar form to the one used in Mack et al. (2002) which treats two kinds of Fourier transforms 
of different normalizations. We here adopt the same normalization to the Fourier transform systematically 
in order to remove uncertainty from numerical calculations, as was pointed out by Lewis (2004) 
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2nW(Koh and Lee, 2000; Mack et al. 2002). In the previous studies the ionization ratio 
in the early Universe was assumed to be a step function of time, i.e. x e = 1 before LSS 
and x e = after that. However, we have to incorporate a correct recombination history in 
order to obtain more accurate theoretical result. We accomplish this by using a numerical 
program RECFAST (Seager et al. 1999) in open use. Thus we take account of the correct 
ionization ratio x e (i]), and we can now rewrite L^(r],\s) = kx e (r])U^ (k). 

Combining Einstein equations with the fluid equations (Ma and Bertschinger 1995, Hu 
and White 1997), we obtain evolution equations of scalar and vector perturbations. We 
evaluated the likelihood functions of WMAP data (Verde et al., 2003) in a wide range of 
parameters of stochastic magnetic field, B\ and %, with other cosmological parameters, 
h,flbh 2 ,fi m h 2 ,n s ,A s , and r in flat Universe models, where VL^h 2 and Q m h 2 are the baryon 
and cold dark matter densities, n s and A s are the spectral index and the amplitude of 
primordial scalar fluctuation, and r is the optical depth. To explore the parameter space, 
we make use of the Markov chain technique (Lewis 2002). We also take account of the SZ 
effect in our analysis. For that, we follow an estimate of Komatsu and Seljak, with o"§ — 0.9 
(Spergel et at. 2003; Komatsu and Seljak 2002). Note that, we consider linear perturbations 
in flat FRW Universe, but we do not consider the effect of gravity wave damping (Caprini 
and Durrer, 2002) which would be relevant for the magnetic fields on super-horizon scales 
if they were generated before the epoch of BBN. Although this might alter the original 
spectrum of primordial magnetic field by the time of LSS in blue power-spectral indices 
— 2 < n (Caprini and Durrer, 2002), one can still allowed to represent the modified spectrum 
in running power-law index. 



3. Result and Discussions 

Our study confirms that the effect of not only the vector mode but the scalar mode 
magnetic field plays the important role in order to solve the potential discrepancy of the 
CMB anisotropies between theory and observation at higher multipoles. 

We show the result of our MCMC analysis with WMAP data in the two parameter plane 
|B A | vs. n B in Fig.l. Although the higher I data from CBI(Mason et al., 2003), ACBAR(Kuo 
et al., 2004), and others may constrain the magnetic field strongly, we did not include those 
data in the present analysis because they have too large error bars. Therefore, the result is 
most sensitive to the WMAP data points at 500< I <900. Shown in Fig.l is the excluded 
region at 2<r(95.4%) C.L as bounded by the thick solid curve, but we could not find the lower 
boundary of the allowed region at the same 2<r(95.4%) C.L. Note that we find a very shallow 
minimum of the reduced x 2 — 1-08 along the thick dashed line which is almost parallel to the 
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upper 2(7 boundary displayed as thick solid line is Fig.l. We can thus obtain the strength 
of primordial magnetic field |B| A <3.9nG(2o") at IMpc for the power-law spectral index 
tib ~ 1.1. This upper limit is the most reliable and newest one to the primordial magnetic 
field because we consider all effects on the CMB anisotropies, i.e. the effect of ionization 
ratio, the SZ effect, and the both scalar and vector mode effects of the magnetic field, for the 
first time, in the present estimate of |B>| and ub- In our numerical estimate of the magnetic 
field parameters, we continued the MCMC analysis until the other cosmological parameters 
converge well to the values listed in table 1. The inferred parameter values are not very 
different from those of Spergel et al. (2002), but the important fact is that we cannot find 
obvious degeneracies of magnetic field parameters with other cosmological parameters. We 
understand this for the following reasons. Primordial magnetic field is constrained by the 
cutoff scale of damping Alfven wave as indicated by Eq.(2) (Jedamzik et al. 1998; Mack 
et al. 2002). Since the cutoff scale is small enough compared with the mulitpoles I <C 
3000 for current available data which we are interested in, the manifestation of this effect 
is remarkably seen as a monotonically increasing contribution from both scalar and vector 
perturbations in the mulutipole region higher then I ~500 as far as |B| A <3.9nG. This 
sensitivity of the CMB power spectrum to the primordial magnetic field differs completely 
from those to the other cosmological parameters, which helps resolve the degeneracies among 
them. 

There is, however, strong degeneracy between magnetic field strength |B A | and power 
spectral index n B . At this moment, we can not resolve this strong degeneracy, because 
the WMAP data are constrained to lower I < 900. The CBI and ACBAR data are in an 
interesting higher multipole region, but the error bars are too larger to resolve the degeneracy. 
Precise data for higher /, where the effect of the primordial magnetic field is especially strong, 
are highly desirable. 

Let us shortly discuss the consistency of our upper limit |B A | < 3.9nG(2<r) with the other 
observational constraint on the magnetic field of the cluster of galaxies. Many astronomical 
observations indicate that the magnetic field strength in the present day cluster of galaxies is 
~ (0.1 — l)fiG. Assuming isotropic collapse for cluster formation, one can straightforwardly 
estimate the magnetic field strength about (1 — 10)nG at the epoch of LSS. If the field 
strength is extremely larger than this critical value ~ (1 — 10)nG, we need some damping 
process, and if it is smaller than ~ (1 — 10)nG, on the other hand, we need amplification 
process. We do not know any viable physical process of damping the magnetic field, but there 
are several amplification processes proposed in literature. We found wide allowed region in 
the two parameter plane (|B^|,n) in Fig.l for the primordial magnetic field which best fits 
the current CMB data. Their parameter values satisfy the observational constraint form the 
cluster of galaxies \H\\ ~ (1 — 10)nG at the epoch of LSS. 
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To summarize, we studied the CMB fluctuation power spectrum by taking account of the 
scalar and vector modes from primordial magnetic field and the SZ effect. The likelihood 
analysis of the WMAP data indicates that the upper limit of magnetic field strength is 
| B A | < 3.9 nG (2a) at 1 Mpc for any power spectral indices % and reasonable values of the 
other cosmological parameters. 
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Fig. 1. — Excluded and allowed regions at 2er(95.4%) C.L. on two parameter plane |B^| vs. 
nB-, where |B>| is the primordial magnetic field strength and ns is the power-law spectral 
index. Thick dashed line is for a shallow minimum of the reduced \ 2 — 1-08, which runs 
through crosses (x) for (|B A |/nG,n B ) = (0. 1,-2.8), (1.9,-1.6), (1.6,-1.3), and (0.076,0.1). 
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Table 1. Calculated ACDM model parameters which best fits the WMAP data with the 
effect of primordial magnetic field taken into account. 



Parameter Mean and 68% C.L. Errors 95% C.L. Errors 



Q b h 2 


o 023+ - 001 


+0.002 
-0.002 


^CDMh 2 


o 1 s+ 002 

u - io -0.01 


+0.03 
-0.03 


n s 


i-ool8:S 


+0.05 
-0.05 


A s 


o.87i8;ll 


+0.20 
-0.20 


T 


o.iil8:8S 


+0.10 
-0.09 


h 


o.7ol8:8 4 


+0.08 
-0.09 



Note. - Calculated corresponding cosmic expansion ages 
are t /Gyr = 13.3+°;1(68% C.L.) and 13.3+°J(95% C.L.). 



